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ABSTRACT: Methane (CH,) wettability of shale is a key parameter which determines pore and reservoir-scale fluid 
distributions, CH, reserves estimation, and ultimate recovery efficiency from shale gas reservoirs. Clay minerals usually fill the 
pore spaces or are adsorbed on the surface of shale rock, thus influencing CH, wettability. However, a systematic investigation 
of the influence of clay on CH, shale wettability is lacking. Herein, we investigated the role of clay, pressure, temperature, and 
salinity on CH, wettability of clay-coated quartz (ie. a well-defined model system for shales). Results indicated that the 
advancing and receding water contact angles for clean, kaolinite-coated, and montmorillonite-coated quartz increased with 
pressure. However, the effect of temperature on wettability is complex, and thus the advancing water contact angle for clean 
quartz increased with temperature while an opposite trend was found for clay-coated quartz. At low temperature (i.e., 300 K), 
clay coating dewetted the quartz surface, while at elevated temperature (i.e., 323 K), clay coating increased the hydrophilicity of 
the quartz surface. Furthermore, kaolinite clay particles demonstrated a stronger influence on quartz wettability than 
montmorillonite particles, both at high and low temperatures. In addition, higher NaCl salinity led to higher advancing water 
contact angles for the aforementioned three solid surfaces. The effect of salinity on CH, wettability is thus intensified in the 
presence of clays. These insights will thus improve the accuracy of CH, reserve estimates and aid methane recovery schemes. 


1. INTRODUCTION 


Shale reserves play an important role in meeting the global 
demands for energy. The significant development in gas 
production from such shale gas reservoirs is due to advances in 
horizontal drilling and hydraulic fracturing technology.'~° 
However, in commercial production, there are challenges in 
evaluating the characteristics of these hydrocarbon reservoirs 
due to the complexity of gas storage and transport processes in 
the micro- and nanoscale pore geometries.” In this context, 
wettability, characterized by the brine contact angles of CH,— 
brine—shale rock (0) is a key parameter that determines CH, 
reserves and production efficiency.” ">? In particular, CH, 
reserves, CH, distribution, and free gas flow in the shale 
pores are strongly affected by the capillary pressure (P.) which 
is a function of contact angle.” 422 For example, for a 
circular capillary tube, capillary pressure can be calculated 
based on the Young—Laplace equation: 


p= 2y cos 0 

Ā r (1) 
where y is the interfacial tension between CH, and brine and r 
is the average pore radius. When 0 < 90°, the shale rock has a 
stronger affinity to brine than CH, and P, could be a driving 
force for CH, recovery as water would imbibe more easily; 
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while if 0 > 90°, the shale rock is more CH,-wet than water- 
wet and P, would tend to resist water imbibition. 

Organic rich shales are mainly composed of consolidated 
clay-sized particles with a high organic content, a complex 
chemical composition that is filled or covered with clay,” 
which complicates contact angle analysis. For example, the clay 
crystal size is much smaller than the usual sessile droplet, 
which prevents standard wettability measurements.”° ~* 
Previous literature reports that kaolinite and montmorillonite 
crystal sizes are smaller than 2 zm, ””° while the usual sessile 
droplet size is around 2000 ym.°'”” To avoid the challenge to 
measure the contact angle of clay using the sessile method, 
Borysenko et al. used a liquid—liquid extraction method to 
validate montmorillonite as a hydrophilic surface and kaolinite 
as a hydrophobic surface.” However, the liquid-liquid 
extraction method could not report the specific contact angle 
value.’ Schmatz et al. reported a more water-wet nature of 
kaolinite based on cryo SEM results.** However, the cryo SEM 
method is destructive.” Ballah et al. measured the wettability 
of clay coating on a glass slide using the sessile method, and all 
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Figure 1. SEM images of clean (a), kaolinite-coated (b, c), and montmorillonite-coated (d, e) quartz used in the experiments. The light color is 


clay particles and the dark color is quartz substrate. 


the contact angles measured in their work were smaller than 
40° (ie, a strong water-wet surface).°° However, all their 
measurements were conducted at ambient conditions, which 
could not represent the real wettability of clay at reservoir 
conditions as pressure and temperature could also influence 
wettability significantly.” A recent study compacted clay 
powders into solid substrates and measured the brine contact 
angles of the CO,/N,/oil—brine—clay systems at reservoir 
conditions.’ They concluded that montmorillonite is oil-wet, 
while kaolinite and illite are water-wet at typical reservoir 
conditions.” However, they did not measure clay wettability in 
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the CH,—brine—clay system, which is related to shale gas 
production.” '* Molecular dynamic simulations, which over- 
come these limitations, predicted that Na-/Ca-montmorillon- 
ite becomes more hydrophobic when exposed to CO, at 
reservoir conditions.” Šolc et al. calculated the contact angle 
of water—air—kaolinite was 105° using molecular dynamic 
simulation." Thus, the inconsistency between experimental 
and simulated results is clear and clay wettability is still open to 
large uncertainty. However, how clay influences the wettability 
of organic rich shales at reservoir conditions has not been 
explored. 
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Figure 2. Effect of pressure and temperature on the (a) advancing and (b) receding brine contact angles of clean/kaolinite-coated/ 
montmorillonite-coated quartz—CH,—brine systems. Note that the aqueous phase is 1.5 wt % NaCl brine. 


In the present study, we thus investigate the CH, wettability 
of clay-coated quartz at reservoir conditions. Note that we use 
the words clay and clay minerals interchangeably in this work. 
To isolate the role of shale constituents, we use quartz, one of 
the most primary components of shale,’’~*’ coated with clay 
particles, instead of real rock samples. Thus, we can elucidate 
the role of clay types on the wettability without perturbations 
from other chemical components in shale. Furthermore, the 
effects of pressure, temperature, and salinity on CH, 
wettability are also presented. 


2. EXPERIMENTAL METHODOLOGY 


To remove any residual contaminants from the quartz surface, 
alpha-quartz substrates (supplied by Dade Quartz Co., Ltd., China, 
with RMS surface roughness of less than 1 nm measured with an JPK 
NanoWizard 4 atomic force microscope (AFM), Figure Sla) were 
carefully cleaned following a standard procedure. First, the substrates 
were washed by acetone and subsequently soaked in piranha solution 
(3 parts of 98 wt % H,SO, and 1 part of aqueous 30 wt % H,O, for 
30 min at 300 K; H,SO, and H,O, were purchased from Sinopharm 
Chemical Reagent Co., Ltd.). Afterward, the immersed substrates 
were covered with aluminum foil and dried in a clean oven for 8 h at 
353 K. A concentration of 2 wt % kaolinite or montmorillonite 
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suspension was prepared respectively through adding kaolinite or 
montmorillonite powders into 1.5 wt % NaCl brine (kaolinite, 
montmorillonite, and sodium chloride were obtained from Sigma- 
Aldrich; NaCl purity was >99 mol %). Note that the chemical 
formula of kaolinite is Al,O,Si,-2H,O and the montmorillonite is a K 
10 catalyst, MgNaAl,(Si,Oj9)3;(OH)¢. A relatively high NaCl salinity 
was used here to impede clay accumulation/swelling and promote 
clay adsorption onto the quartz surface.°°*' The total salinities of clay 
suspensions in DI water and in 1.5 wt % NaCl brine were measured 
using Starter3100 M at room temperature (see Table S1). It was 
found that the total salinity of clay suspensions is extremely low, 
showing no apparent difference with the salinity of pure DI water. A 
concern may be that such a high Na* concentration would exchange 
with in situ divalent cations from clay minerals. However, as we show 
in Tables S1 and S2, the ion salinities from in situ clay were extremely 
low and the wettabilities for Na* and Mg”* were very similar. Thus, it 
was assumed that the ion exchange would not influence wettability in 
this work. The particle size distributions (PSD) of the kaolinite and 
montmorillonite were measured thrice via a Zetasizer Nano 
instrument (Malvern, U.K.). The average particle sizes for kaolinite 
and montmorillonite were 1050 and 900 nm, respectively. The total 
organic content (TOC) for kaolinite and montmorillonite were 335 
and 3093 mg/kg, respectively, measured by an Elab-TOC instrument. 
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Figure 3. Advancing brine contact angles of CH,—clean/kaolinite-coated/montmorillonite-coated quartz—brine systems at 10 MPa and 323 K. 


The kaolinite and montmorillonite suspensions were stirred for 2 h 
at 323 K. Afterward, the clean quartz substrates were aged in the 2 wt 
% kaolinite or montmorillonite suspensions for 6 h, and then the aged 
substrates were covered with clean aluminum foil and dried for 6 h. 
The total 12 h duration would ensure a stable coating. The 
morphological properties of clean and coated substrates were 
characterized using AFM (note that RMS roughness values of clean 
and kaolinite-coated quartz were measured by a JPK NanoWizard 4 
atomic force microscope, while RMS of montmorillonite quartz was 
measured via a Bruker Multi Mode 8 atomic force microscope, please 
see Figure S1) and SEM measurements by a Philips XL 30 SEM 
(Figure 1). The RMS roughness values for clean, kaolinite-coated, and 
montmorillonite-coated quartz were 483.7 pm, 629.8 nm, and 863 
nm, respectively. However, at this scale, the effect of roughness on 
contact angle is insignificant.°* Based on the analysis of SEM images 
using MATLAB, the average kaolinite and montmorillonite coverages 
were 37.48% and 39.02%, respectively (see Figure $2). Furthermore, 
the SEM images demonstrated a more compact coating on the 
montmorillonite-coated surface. The advancing (6,) and receding (0,) 
brine contact angles of CH,—brine—clean/clay-coated quartz systems 
were then measured by a Kriiss DSA 100 instrument using the tilted 
plate method? as described previously in detail”? All the contact 
angle measurements were finished within the initial 30 s. We did not 
observe the contact angle for a long time in this work, because the 
droplet could destabilize and displace the clay coating after long time 
contact. However, a future work would investigate long-term clay 
wettability using a recent developed LbL coating method.** 
Furthermore, we did not consider the influence from CH,—brine 
dissolution because of a very low dissolution capacity between brine 
and CH,.°° Note that in shale gas extraction, 0, corresponds to CH, 
recovery, while 0, is related to initial CH,/water distribution.°° 


3. RESULTS AND DISCUSSION 


3.1. Effect of Pressure and Temperature on Contact 
Angles. The effects of pressure and temperature on the brine 
contact angles were measured for the various substrates at 
multiple pressures (0.1 to 20 MPa) and temperatures (300 and 
323 K), see Figure 2. 

Both 6, and 0, for all substrates clearly increased with 
pressure, consistent with literature data measured for several 
different gases on various substrates.7°39435%57-5? Ror 
example, at 300 K, when pressure increased from 0.1 to 20 
MPa, @, for clean, kaolinite-coated, and montmorillonite- 
coated quartz increased by 30°, 29°, and 21°, respectively, 
while 0, for the above three surfaces increased by 22°, 25°, and 
15°, respectively. 
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Noteworthy was the effect of temperature on the wettability. 
0, for clean quartz increased with temperature, consistent with 
literature data reported for quartz,” hydrophilic dolomite,*’ 
and shale,” while 0, for clay-coated quartz decreased with 
temperature, consistent with literature data reported for 
mica,” hydrophobic dolomite,’ and 3 wt % TOC shale.” 
Specifically, at 15 MPa, when temperature increased from 300 
to 323 K, @, changes for clean, kaolinite-coated, and 
montmorillonite-coated quartz were 9°, —13°, and —8°, 
respectively; at 20 MPa, for the same temperature increase, 
0, changes for the above three surfaces were 5°, —24°, and 
—6°, respectively. 

Physically, the influence of pressure on wettability can be 
attributed to the higher nonaqueous phase density, >” the 
consequently stronger intermolecular interactions between 
CH, and the solid surface due to the increased molecular 
density, °°" and the corresponding decrease in the 
interfacial tension between the solid surface and CH,.* The 
underlying mechanism of the temperature influence on 
wettability can be explained by the sharp-kink approxima- 


tion, >” as shown in eq 2 below. 
ApI 
cos 6 = —1 + 2 
Ng (2) 


where Ap is the approximate density difference between the 
adsorbed brine film and CH4 y is the interfacial tension 
between brine and CH, and I is the van der Waals potential 
integral. 

Clearly, both yj, and Ap increased with temperature.” 
However, the relative increase in Ap with temperature was 
more significant than that of ie which reduces the contact 
angle; the van der Waals potential decreases with temper- 
ature,“ which increases the contact angle. Thus, in both 
scenarios, a 0 increase or decrease with temperature is possible. 
For clean quartz, the temperature-driven influence on van der 
Waals interaction was more effective, resulting in a contact 
angle increase with temperature.’ For clay-coated quartz, this 
van der Waals interaction might be relatively weak, and thus 
the contact angle decreased with temperature. 

3.2. Effect of Brine Salinity on Contact Angles. All 0, 
increased with salinity as depicted in Figure 3. Specifically, 
when NaCl salinity increased from 1.5 to 15 wt %, 0, for clean, 
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Table 1. Zeta Potential and pH Data for Various Systems at Two Different Temperatures in This Work 


sample T (K) mean zeta potential (mV) standard deviation (mV) mean pH standard deviation 
1.5 wt % NaCl 300 —0.86 1.01 6.45 
1.5 wt % NaCl 323 1.21 0.786 
S wt % NaCl 300 1.91 4.94 6.25 0.03 
S wt % NaCl 323 3.49 6.65 
2 wt % kaolinite + 1.5 wt % NaCl 300 -14.7 0.436 6.25 0.58 
2 wt % montmorillonite + 1.5 wt % NaCl 300 21.9 0.574 4.74 0.23 
2 wt % kaolinite + 1.5 wt % NaCl 323 -12.6 0.872 
2 wt % montmorillonite + 1.5 wt % NaCl 323 -19.6 0.967 
2 wt % kaolinite + 5 wt % NaCl 300 -7.51 7.88 5.61 0.02 
2 wt % montmorillonite + 5 wt % NaCl 300 -10.1 6.3 3.48 0.01 
2 wt % kaolinite + 5 wt % NaCl 323 —3.04 6.58 
2 wt % montmorillonite + 5 wt % NaCl 323 —2.14 3.59 
2 wt % kaolinite + 10 wt % NaCl 323 2.37 9.31 
2 wt % montmorillonite + 10 wt % NaCl 323 3.48 8.97 


kaolinite-coated, and montmorillonite-coated quartz increased 
by 6°, 15°, and 10°, respectively. 

This result is consistent with literature data for NaCl—oil— 
sandstone,” NaCl—glass/kaolinite-coated glass—crude oil, 
NaCl—oil—kaolinite-coated glass,” NaCl—CO,—mica,*8>”"”° 
and NaCl—CO,—quartz**” systems. The effect of salinity on 
the advancing water contact angles can be explained with the 
Derjaguin—Landau—Verwey—Overbeek (DLVO) theory as 
discussed in our previous work. In order to shed more 
light on the underlying mechanism of salinity influence on the 
contact angles, we also measured the corresponding zeta 
potentials (using a Zetasizer Nano instrument, Malvern, U.K.) 
and pH values (using a RHA-SHKY pH meter model 8601 
instrument), see Table 1. Every zeta-potential measurement 
was repeated 18 times, and every pH measurement was 
repeated 3 times. 

Higher temperatures led to less negative zeta potentials, 
while the addition of kaolinite or montmorillonite reduced the 
zeta potential. Furthermore, clay suspensions are not stable 
and uniform solutions, and this intrinsic property could cause a 
significant variation in the measured properties.’' 

In addition, montmorillonite suspensions were more acidic 
than the corresponding kaolinite suspensions, while a higher 
NaC] salinity led to a lower pH value for both montmorillonite 
and kaolinite suspensions. Higher pH values led to more 
negative zeta potentials” and thus lower contact angles. 
However, it could not be explained why the contact angle for 2 
wt % montmorillonite + 1.5 (or 5) wt % NaCl brine sample 
was smaller than that for 2 wt % kaolinite + 1.5 (or 5) wt % 
NaCl brine sample at 300 K just based on pH. We hypothesize 
that the surface microstructural evolution and (clay) hydration 
of the coated surface also played an important role, causing the 
different wettability behaviors for kaolinite- and montmor- 
illonite-coated surfaces, besides the solution property.°””* 
Clearly, the contact angle strongly correlated with the zeta 
potential, as shown in Table 2 and Figure 4. 

At 323 K, a higher salinity led to a larger zeta potential, thus 
higher contact angles. The zeta potential for kaolinite was 
smaller than that for montmorillonite at the same NaCl salinity 
(S and 10 wt %, respectively), and thus the contact angle for 
kaolinite was smaller than montmorillonite. However, note 
that the zeta potential was measured at ambient conditions, 
while the contact angle was measured at high pressure, 
although a correlation has been proposed previously.’* 
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Table 2. Used Data for Plotting Figure 4 


zeta potential advancing contact 
solution system mV) angle (deg) 
2 wt % montmorillonite + —19.6 15 
1.5 wt % NaCl 
2 wt % kaolinite + 1.5 wt % NaCl —12.6 15 
2 wt % kaolinite + 5 wt % NaCl —3.04 22 
2 wt % montmorillonite + —2.14 23 
5 wt % NaCl 
2 wt % kaolinite + 10 wt % NaCl 2.37 24 
2 wt % montmorillonite + 3.48 26 
10 wt % NaCl 
30 
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Figure 4. Relationship between advancing contact angle and the 
corresponding zeta potential at 10 MPa, 323 K, and various NaCl 
salinities. 


In summary, the increasing salinity increases the ionic 
strength, which better screens the surface charge.°! Thus, the 
surface potential approaches zero and, in some circumstances, 
positive values. As a result, surface polarity and water affinity 
decreases, and hence, an increasing water contact angle results. 

3.3. Effect of Clay on Quartz Wettability. Based on the 
results and discussion in Sections 3.land 3.2, it is easy to 
conclude that the effect of clay on quartz wettability is 
extremely complex. For example, the effect of clay coating on 
quartz wettability was different at low and high temperatures, 
see Figure 2. At low temperature (i.e, 300 K), kaolinite or 
montmorillonite coating dewetted the quartz, which could be 
attributed to the lowered pH and intensified acidic properties, 
thus less negative zeta potential,” and a larger contact angle,”* 
while at high temperature (ie, 323 K), these coatings 
intensified quartz hydrophilicity. 
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Furthermore, kaolinite always had a stronger influence on 
quartz wettability than montmorillonite. For example, at 20 
MPa and 300 K, 0, values for clean, kaolinite-coated, and 
montmorillonite-coated quartz were 30°, 47°, and 34°, 
respectively, while at 20 MPa and 323 K, 0, values for the 
above three surfaces were 35°, 23°, and 28°. The low 
temperature results (300 K) are consistent with literature 
reports that (a) kaolinite and montmorillonite clays were more 
oil-wet than quartz at 298 K and 0.1 MPa,” that (b) quartz 
had a stronger water adsorption capacity than kaolinite at 298 
K,” and that (c) brine contact angles of oil—-water—sandstone 
systems increased with clay content on the sandstone surface at 
ambient conditions (i.e., 0.1 MPa and 287 Ki? Furthermore, 
the results for the clay influence at 323 K are consistent with 
literature reports where kaolinite reduced water repellency of 
silica sand at 323 K more efficiently than montmorillonite.”° 
Although a recent study reported clay wettability for CO,/N,/ 
oil—brine—clay systems at reservoir conditions,” we did not 
find any literature related to clay influence on shale wettability 
at high temperature and high pressure, noting that reservoir 
conditions (i.e., elevated pressure and temperature) are critical 
for gas behavior at the pore-scale,” which is again vital for the 
evaluation of CH, reserves and recovery. 

In addition, the effect of salinity on CH, wettability was 
intensified in the presence of clays at 10 MPa and 323 K, see 
Figure 3. Specifically, when the NaCl salinity increased from 
1.5 to 15 wt % for clean quartz, 0, increased by 6° from 25° to 
31°; compared with this, salinity had a stronger influence on 0, 
for the kaolinite-coated (e.g., 0, increased by 10° from 15° to 
25°) and montmorillonite-coated (e.g, 0, increased by 15° 
from 15° to 30°) surfaces. Mechanistically, clay particles have 
negatively charged surfaces at typical reservoir conditions, ® 
consistent with our zeta-potential results; the isoelectric points 
for quartz, kaolinite, and montmorillonite are pH = 3,” 6, and 
6.7.” Montmorillonite has a greater number of structurally 
charged sites than kaolinite,’” which results in a stronger 
attraction to positive charges (e.g, Na*). Thus, the surface 
potential for montmorillonite-coated quartz approaches zero 
more quickly. Consequently, 9, for montmorillonite-coated 
quartz increased more significantly than 0, for both kaolinite- 
coated quartz and clean quartz. 


4. CONCLUSIONS 


CH,-—shale wettability is a key parameter for predicting CH, 
reserves and distribution.” ~"? However, due to the complexity 
of the shale structure and composition, CH,—shale wettability 
evaluation is still open to uncertainties. Specifically, clay 
minerals, as one of the abundant constituents of shale 
formations, significantly increase the complexity of CH,— 
shale wettability. Furthermore, a systematic investigation of 
clay effects on CHy,-shale wettability in a representative 
system is still scarce. To elucidate how clay influences shale 
wettability, we used a well-defined model (i.e, clay-coated 
quartz) to investigate the role of clay, pressure, temperature, 
and salinity on the CH, wettability. The following conclusions 
were reached: 

(1) The quartz substrates became more water-wet at high 
temperature (i.e, 323 K) when coated by clays, while they 
became more CH,-wet at low temperatures (i.e, 300 K). 
Kaolinite had a stronger influence on quartz wettability than 
montmorillonite with respect to temperature variations. We 
did not find any previous reports to compare the influence of 


793 


clay coating or adsorption on shale and quartz wettability at 
different temperatures, especially at reservoir pressures. 

(2) Elevated pressures increased the hydrophobicity of all 
three tested surfaces (i.e clean, kaolinite-coated, and 
montmorillonite-coated quartz). 

(3) Elevated temperatures dewetted the clean quartz, while 
they increased the hydrophilicity of clay-coated quartz. 

(4) High NaCl salinity reduced surface hydrophilicity. This 
salinity effect was intensified in the presence of clays. 
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